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A measurement of the cross section for the inclusive production of isolated prompt photons in 
pp collisions at a centre-of-mass energy y^s — 7 TeV is presented. The measurement covers the 
pseudorapidity ranges I?;'''! < 1.37 and 1.52 < < 1.81 in the transverse energy range 15 < < 
100 GeV. The results are based on an integrated luminosity of 880 nb~^, collected with the ATLAS 
detector at the Large Hadron Collider. Photon candidates are identified by combining information 
from the calorimeters and from the inner tracker. Residual background in the selected sample is 
estimated from data based on the observed distribution of the transverse isolation energy in a narrow 
cone around the photon candidate. The results are compared to predictions from next-to-leading 
order perturbative QCD calculations. 
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I. INTRODUCTION 

Prompt photon production at hadron colliders pro- 
vides a handle for testing perturbative QCD (pQCD) 
predictions [1, 2]. Photons provide a colorless probe of 
quarks in the hard partonic interaction and the subse- 
quent parton shower. Their production is directly sen- 
sitive to the gluon content of the proton through the 
qg — > 57 process, which dominates at leading order 
(LO). The measurement of the prompt photon produc- 
tion cross section can thus be exploited to constrain the 
giuon density function [3, 4]. Furthermore, photon identi- 
fication is important for many physics signatures, includ- 
ing searches for Higgs boson [5-7] , graviton decays [8] to 
photon pairs, decays of excited fermions [9], and decays 
of pairs of supersymmetric particles characterized by the 
production of two energetic photons and large missing 
transverse energy [10-12]. 

Prompt photons include both "direct" photons, which 
take part in the hard scattering subprocess (mostly 
quark-gluon Compton scattering, qg — > 57, or quark- 
antiquark annihilation, qq — > 97), and "fragmenta- 
tion" photons, which arc the result of the fragmenta- 
tion of a high-pT parton [13, 14]. In this analysis, 
an isolation criterion is applied based on the amount 
of transverse energy inside a cone of radius R = 

\J {rj — 11^)^ + (0 — (lf<)^ ~ 0.4 centered around the pho- 
ton direction in the pseudorapidity (77) and azimuthal 
angle (0) plane [15]. After the isolation requirement is 
applied the relative contribution to the total cross section 
from fragmentation photons decreases, though it remains 
non- negligible especially at low transverse energies [14]. 
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The isolation requirement also significantly reduces the 
main background of non-prompt photon candidates from 
decays of energetic tt*^ and r/ mesons inside jets. 

Early studies of prompt photon production were car- 
ried out at the ISR collider [16, 17]. Subsequent stud- 
ies, for example [18-20], further established prompt pho- 
tons as a useful probe of parton interactions. More re- 
cent measurements at hadron colliders were performed at 
the Tcvatron, in pp collisions at a centre-of-mass energy 
y/s ~ 1.96 TeV. The measurement by the DO Collabo- 
ration [21] is based on 326 pb~^ and covers a pseudo- 
rapidity range \rf\ < 0.9 and a transverse energy range 
23 < < 300 GeV, while the measurement by the 
CDF Collaboration [22] is based on 2.5 fb~^ and covers a 
pseudorapidity range |77'''| < 1.0 and a transverse energy 
range 30 < < 400 GeV. Both DO and CDF measure 
an isolated prompt photon cross section in agreement 
with next-to-leading order (NLO) pQCD calculations, 
with a slight excess seen in the CDF data between 30 
and 50 GcV. Measurements of the inclusive prompt pho- 
ton production cross section have also been performed 
in ep collisions, both in photoproduction and deep in- 
elastic scattering, by the HI [23, 24] and ZEUS [25, 26] 
Collaborations. The most recent measurement of the in- 
clusive isolated prompt photon production was done with 
2.9 pb-i at 7 TcV by the CMS Collaboration [27]. 

That measurement, which covers 21 < < 300 GeV 
and [77'''] < 1.45, is in good agreement with NLO predic- 
tions for the full E^ range. 

This paper describes the extraction of a signal of iso- 
lated prompt photons using 880 nb~^ of data collected 
with the ATLAS detector at the Large Hadron Collider 
(LHC). A measurement of the production cross section 
in pp collisions at ^/s = 7 TeV is presented, in the pseu- 
dorapidity ranges Ir]^] < 0.6, 0.6 < < 1.37 and 
1.52 < \ri'^\ < 1.81, for photons with transverse energies 
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between 15 GeV and 100 GeV. 

The paper is organized as follows. The detector is de- 
scribed in Section II, followed by a summary of the data 
and the simulated samples used in the measurement in 
Section 111. Section IV introduces the theoretical predic- 
tions to which the measurement is compared. Section V 
describes the photon reconstruction and identification al- 
gorithms; their performance is given in Section VI. Sec- 
tion VII describes the methods used to estimate the resid- 
ual background in the data and to extract the prompt 
photon signal, followed by a discussion of the data correc- 
tions for the cross section measurement in Section VIII. 
The sources of systematic uncertainties on the cross sec- 
tion measurement are discussed in Section IX. Section X 
contains the main experimental results and the compar- 
ison of the observed cross sections with the theoretical 
predictions, followed by the conclusions in Section XI. 



II. THE ATLAS DETECTOR 

The ATLAS detector is described in detail in Refs. [28] 
and [29]. For the measurement presented in this paper, 
the calorimeter, with mainly its electromagnetic section, 
and the inner detector arc of particular relevance. 

The inner detector consists of three subsystems: at 
small radial distance r from the beam axis (50.5 < 
r < 150 mm), pixel silicon detectors are arranged in 
three cylindrical layers in the barrel and in three disks 
in each end-cap; at intermediate radii (299 < r < 560 
mm), double layers of single-sided silicon microstrip de- 
tectors are used, organized in four cylindrical layers in 
the barrel and nine disks in each end-cap; at larger radii 
(563 < r < 1066 mm), a straw tracker with transition 
radiation detection capabilities divided into one barrel 
section (with 73 layers of straws parallel to the beam 
line) and two end-caps (with 160 layers each of straws 
radial to the beam line) is used. These three systems are 
immersed in a 2 T axial magnetic field provided by a su- 
perconducting solenoid. The inner detector has full cov- 
erage in (j). The silicon pixel and microstrip subsystems 
cover the pseudorapidity range |?7| < 2.5, while the transi- 
tion radiation tracker (TRT) acceptance is limited to the 
range \r]\ < 2.0. The inner detector allows an accurate 
reconstruction of tracks from the primary proton-proton 
collision region, and also identifies tracks from secondary 
vertices, permitting the efficient reconstruction of pho- 
ton conversions in the inner detector up to a radius of 
sa 80 cm. 

The electromagnetic calorimeter is a lead-liquid argon 
(Pb-LAr) sampling calorimeter with an accordion geom- 
etry. It is divided into a barrel section, covering the pseu- 
dorapidity region \r]\ < 1.475, and two end-cap sections, 
covering the pseudorapidity regions 1.375 < < 3.2. 
It consists of three longitudinal layers. The first one, 
with a thickness between 3 and 5 radiation lengths, is 
segmented into high granularity strips in the r] direction 
(width between 0.003 and 0.006 depending on ij, with the 



exception of the regions 1.4 < \r]\ < 1.5 and |?7| > 2.4), 
sufficient to provide an event-by-event discrimination be- 
tween single photon showers and two overlapping showers 
coming from a tt" decay. The second layer of the electro- 
magnetic calorimeter, which collects most of the energy 
deposited in the calorimeter by the photon shower, has 
a thickness around 17 radiation lengths and a granular- 
ity of 0.025 X 0.025 in 77 x (corresponding to one cell). 
A third layer, with thickness varying between 4 and 15 
radiation lengths, is used to correct leakage beyond the 
calorimeter for high energy showers. In front of the ac- 
cordion calorimeter a thin presampler layer, covering the 
pseudorapidity interval |?7| < 1.8, is used to correct for 
energy loss before the calorimeter. The sampling term a 
of the energy resolution {a{E)/ E k, a/ E [GeV]) varies 
between 10% and 17% as a function of |?7| and is the 
largest contribution to the resolution up to about 200 
GeV, where the global constant term, estimated to be 
0.7% [30], starts to dominate. 

The total amount of material before the first active 
layer of the electromagnetic calorimeter (including the 
presampler) varies between 2.5 and 6 radiation lengths 
as a function of pseudorapidity, excluding the transition 
region (1.37 < \'q\ < 1.52) between the barrel and the 
end-caps, where the material thickness increases to 11.5 
radiation lengths. The central region {\ri\ < 0.6) has 
significantly less material than the outer barrel (0.6 < 
1 77! < 1.37), which motivates the division of the barrel 
into two separate regions in pseudorapidity. 

A hadronic sampling calorimeter is located beyond the 
electromagnetic calorimeter. It is made of steel and scin- 
tillating tiles in the barrel section (|?7| < 1.7), with depth 
around 7.4 interaction lengths, and of two wheels of cop- 
per and liquid argon in each end-cap, with depth around 
9 interaction lengths. 

A three-level trigger system is used to select events con- 
taining photon candidates during data taking. The first 
level trigger (level- 1) is hardware based: using a coarser 
cell granularity (0.1 x 0.1 in 77 x 0) than that of the elec- 
tromagnetic calorimeter, it searches for electromagnetic 
clusters within a fixed window of size 0.2x0.2 and retains 
only those whose total transverse energy in two adjacent 
cells is above a programmable threshold. The second and 
third level triggers (collectively referred to as the "high- 
level" trigger) are implemented in software. The high- 
level trigger exploits the full granularity and precision 
of the calorimeter to refine the level- 1 trigger selection, 
based on improved energy resolution and detailed infor- 
mation on energy deposition in the calorimeter cells. 



III. COLLISION DATA AND SIMULATED 
SAMPLES 

A. Collision Data 

The measurement presented here is based on proton- 
proton collision data collected at a centre-of-mass energy 
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y/s = 7 TeV between April and August 2010. Events 
in which the calorimeters or the inner detector are not 
fully operational, or show data quality problems, are ex- 
cluded. Events are triggered using a single-photon high- 
level trigger with a nominal transverse energy thresh- 
old of 10 GeV, seeded by a level-1 trigger with nomi- 
nal threshold equal to 5 GeV. The selection criteria ap- 
plied by the trigger on shower shape variables computed 
from the energy profiles of the showers in the calorime- 
ters are looser than the photon identification criteria ap- 
plied in the offline analysis, and allow ATLAS to reach 
a plateau of constant efficiency close to 100% for true 
prompt photons with > 15 GeV and pseudorapid- 
ity |?7''| < 1.81. In addition, samples of minimum-bias 
events, triggered by using two sets of scintillator counters 
located at z = ±3.5 m from the collision centre, are used 
to estimate the single-photon trigger efficiency. The total 
integrated luminosity of the sample passing data quality 
and trigger requirements amounts to (880 ± 100) nb^^. 

In order to reduce non-collision backgrounds, events 
are required to have at least one reconstructed primary 
vertex consistent with the average beam spot position 
and with at least three associated tracks. The efficiency 
of this requirement is expected to be greater than 99.9% 
in events containing a prompt photon with Ej, > 15 GeV 
and lying within the calorimeter acceptance. The total 
number of selected events in data after this requirement 
is 9.6 million. The remaining amount of non-collision 
background is estimated using control samples collected 
with dedicated, low threshold triggers that are activated 
in events where either no proton bunch or only one of the 
two beams crosses the interaction region. The estimated 
contribution to the final photon sample is less than 0.1% 
and is therefore neglected. 



B. Simulated events 

To study the characteristics of signal and back- 
ground events, Monte Carlo (MC) samples are gener- 
ated using PYTHIA 6.4.21 [31], a leading-order parton- 
shower MC generator, with the modified leading order 
MRST2007 [32] parton distribution functions (PDFs). It 
accounts for QED radiation emitted off quarks in the 
initial state (ISR) and in the final state (FSR). PYTHIA 
simulates the underlying event using the multiple-parton 
interaction model, and uses the Lund string model for 
hadronisation [33]. The event generator parameters are 
set according to the ATLAS MC09 tune [34], and the 
detector response is simulated using the GEANT4 pro- 
gram [35]. These samples are then reconstructed with 
the same algorithms used for data. More details on the 
event generation and simulation infrastructure are pro- 
vided in Ref. [36]. For the study of systematic uncer- 
tainties related to the choice of the event generator and 
the parton shower model, alternative samples are also 
generated with HERWIG 6.5 [37]. This generator also uses 
LO pQCD matrix elements, but has a different parton 



shower model (angle-ordered instead of px-ordered), a 
different hadronisation model (the cluster model) and a 
different underlying event model, which is generated us- 
ing the JIMMY package [38] with multiple parton interac- 
tions enabled. The HERWIG event generation parameters 
are also set according to the MC09 tune. 

To study the main background processes, simulated 
samples of all relevant 2— )-2 QCD hard subprocesses in- 
volving only partons are used. The prompt photon con- 
tribution arising from initial and final state radiation 
emitted off quarks is removed from these samples in stud- 
ies of the background. 

Two different simulated samples are employed to study 
the properties of the prompt photon signal. The first 
sample consists of leading-order 7-jet events, and con- 
tains primarily direct photons produced in the hard sub- 
processes qg qj and qq 177. The second signal sam- 
ple includes ISR and FSR photons emitted off quarks in 
all 2— >-2 QCD processes involving only quarks and glu- 
ons in the hard scatter. This sample is used to study 
the contribution to the prompt photon signal by photons 
from fragmentation, or from radiative corrections to the 
direct process, that are less isolated than those from the 
LO direct processes. 

Finally, a sample oi W ev simulated events is used 
for the efficiency and purity studies involving electrons 
from W decays. 



IV. THEORETICAL PREDICTIONS 

The expected isolated prompt photon production cross 
section as a function of the photon transverse energy E^ 
is calculated with the JETPHOX Monte Carlo program [13], 
which implements a full NLO QCD calculation of both 
the direct and the fragmentation contributions to the to- 
tal cross section. In the calculation performed for this 
measurement, the total transverse energy carried by the 
partons inside a cone of radius i? = 0.4 in the 77 — space 
around the photon direction is required to be less than 
4 GeV. The NLO photon fragmentation function [39] and 
the CTEQ 6.6 parton density functions [40] provided by 
the LHAPDF package [41] are used. The nominal renor- 
malization {(J-r), factorization (fip) and fragmentation 
(fif) scales are set to the photon transverse energy E^. 
Varying the CTEQ PDFs within the 68% C.L. intervals 
causes the cross section to vary between 5% and 2% as 
Et increases between 15 and 100 GeV. The variation of 
the three scales independently between 0.5 and 2.0 times 
the nominal scale changes the predicted cross section by 
20% at low Et and 10% at high E^, while the varia- 
tion of the isolation requirement between 2 and 6 GeV 
changes the predicted cross section by no more than 2%. 
The MSTW 2008 PDFs [42] are used as a cross-check of 
the choice of PDF. The central values obtained with the 
MSTW 2008 PDFs are between 3 and 5% higher than 
those predicted using the CTEQ 6.6 PDFs. 

The NLO calculation provided by JETPHOX predicts a 
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cross section at parton level, which does not include ef- 
fects of hadronisation nor the underlying event and pileup 
(i.e. multiple proton-proton interactions in the same 
bunch crossing). The non-perturbative effects on the 
cross section due to hadronisation are evaluated using the 
simulated PYTHIA and HERWIG signal samples described in 
Section IIIB, by evaluating the ratio of the cross section 
before and after hadronisation and underlying event sim- 
ulation. The ratios are estimated to be within 1% (2%) 
of unity in PYTHIA (HERWIG) for all and 77 regions 
under study. To account for the effect of the underly- 
ing event and pileup on the measured isolation energy, a 
correction to the photon transverse isolation energy mea- 
sured in data is applied, using a procedure described in 
Section VC. 



V. PHOTON RECONSTRUCTION, 
IDENTIFICATION AND ISOLATION 

A. Photon reconstruction and preselection 

Photon reconstruction is seeded by clusters in the elec- 
tromagnetic calorimeter with transverse energies exceed- 
ing 2.5 GeV, measured in projective towers of 3x5 cells 
in 7/ X in the second layer of the calorimeter. An at- 
tempt is made to match these clusters with tracks that 
arc reconstructed in the inner detector and extrapolated 
to the calorimeter. Clusters without matching tracks are 
directly classified as "unconverted" photon candidates. 
Clusters with matched tracks arc considered as elec- 
tron candidates. To recover photon conversions, clusters 
matched to pairs of tracks originating from reconstructed 
conversion vertices in the inner detector are considered 
as "converted" photon candidates. To increase the recon- 
struction efficiency of converted photons, conversion can- 
didates where only one of the two tracks is reconstructed 
(and does not have any hit in the innermost layer of the 
pixel detector) are also retained [29, 30]. 

The final energy measurement, for both converted and 
unconverted photons, is made using only the calorimeter, 
with a cluster size that depends on the photon classifica- 
tion. In the barrel, a cluster corresponding to 3x 5 (?7 x 0) 
cells in the second layer is used for unconverted photons, 
while a cluster of 3 x 7 (77 x 0) cells is used for converted 
photon candidates (to compensate for the opening be- 
tween the conversion products in the 4> direction due to 
the magnetic field). In the end-cap, a cluster size of 5x5 
is used for all candidates. A dedicated energy calibra- 
tion [29] is then applied separately for converted and 
unconverted photon candidates to account for upstream 
energy loss and both lateral and longitudinal leakage. 

Photon candidates with calibrated transverse energies 
(E^) above 15 GeV are retained for the successive anal- 
ysis steps. To minimise the systematic uncertainties re- 
lated to the efficiency measurement at this early stage of 
the experiment, the cluster barycenter in the second layer 
of the electromagnetic calorimeter is required to lie in the 



pseudorapidity region [77 '''j < 1.37, or 1.52 < |77'>'| < 1.81. 
Photon candidates with clusters containing cells over- 
lapping with few problematic regions of the calorimeter 
readout are removed. After the above preselection, 1.3 
million photon candidates remain in the data sample. 



B. Photon identification 

Shape variables computed from the lateral and longitu- 
dinal energy profiles of the shower in the calorimeters are 
used to discriminate signal from background. The exact 
definitions of the discriminating variables are provided 
in Appendix A. Two sets of selection criteria (denoted 
"loose" and "tight" ) are defined, each based on indepen- 
dent requirements on several shape variables. The se- 
lection criteria do not depend on the photon candidate 
transverse energy, but vary as a function of the photon 
reconstructed pseudorapidity, to take into account varia- 
tions in the total thickness of the upstream material and 
in the calorimeter geometry. 



1. Loose identification criteria 

A set of loose identification criteria for photons is de- 
fined based on independent requirements on three quan- 
tities: 

• the leakage i?had in the first layer of the hadronic 
compartment beyond the electromagnetic cluster, 
defined as the ratio between the transverse energy 
deposited in the first layer of the hadronic calorime- 
ter and the transverse energy of the photon candi- 
date; 

• the ratio i?,, between the energy deposits in 3 x 7 
and 7x7 cells in the second layer of the electro- 
magnetic calorimeter; 

• the RMS width W2 of the energy distribution 
along 77 in the second layer of the electromagnetic 
calorimeter. 

True prompt photons are expected to have small hadronic 
leakage (typically below 1-2%) and a narrower energy 
profile in the electromagnetic calorimeter, more concen- 
trated in the core of the cluster, with respect to back- 
ground photon candidates from jets. 

The loose identification criteria on -Rhad, Rri and W2 
are identical for converted and unconverted candidates. 
They have been chosen, using simulated prompt pho- 
ton events, in order to obtain a prompt photon effi- 
ciency, with respect to reconstruction, rising from 97% 
at e:^ = 20 GeV to above 99% for e:^ > 40 GeV for 
both converted and unconverted photons [30] . The num- 
ber of photon candidates in data passing the preselection 
and loose photon identification criteria is 0.8 million. 
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2. Tight identification criteria 

To further reject the background, the selection require- 
ments on the quantities used in the loose identification 
are tightened. In addition, the transverse shape along 
the (j) direction in the second layer (the variable R^, com- 
puted from the ratio between the energy deposits in 3 x 3 
and 3x7 cells) and the shower shapes in the first layer of 
the calorimeter are examined. Several variables that dis- 
criminate single photon showers from overlapping nearby 
showers (in particular those which originate from neutral 
meson decays to photon pairs) are computed from the 
energy deposited in the first layer: 

• the total RMS width w^.tot of the energy distribu- 
tion along Tj; 

• the asymmetry -Eratio between the first and second 
maxima in the energy profile along rj; 

• the energy difference AE between the second max- 
imum and the minimum between the two maxima; 

• the fraction inside of the energy in seven strips cen- 
tered (in rj) around the first maximum that is not 
contained in the three core strips; 

• the RMS width Ws.3 of the energy distribution com- 
puted with the three core strips. 

The first variable rejects candidates with wide show- 
ers consistent with jets. The second and third variables 
provide rejection against cases where two showers give 
separated energy maxima in the first layer. The last two 
variables provide rejection against cases where two show- 
ers are merged in a wider maximum. 

The tight selection criteria are optimised indepen- 
dently for unconverted and converted photons to account 
for the quite different developments of the showers in each 
case. They have been determined using samples of simu- 
lated signal and background events prior to data taking, 
aiming to obtain a average efficiency of 85% with respect 
to reconstruction for true prompt photons with trans- 
verse energies greater than 20 GeV [30]. About 0.2 mil- 
lion photon candidates are retained in the data sample 
after applying the tight identification requirements. 

C. Photon transverse isolation energy 

An experimental isolation requirement, based on the 
transverse energy deposited in the calorimeters in a cone 
around the photon candidate, is used in this measure- 
ment to identify isolated prompt photons and to further 
suppress the main background from tt*^ (or other neu- 
tral hadrons decaying in two photons), where the 7r° is 
unlikely to carry the full original jet energy. The trans- 
verse isolation energy {E}^°) is computed using calorime- 
ter cells from both the electromagnetic and hadronic 
calorimeters, in a cone of radius 0.4 in the rj ~ (j) space 



around the photon candidate. The contributions from 
5x7 electromagnetic calorimeter cells in the rj — (j) space 
around the photon barycenter are not included in the 
sum. The mean value of the small leakage of the photon 
energy outside this region, evaluated as a function of the 
photon transverse energy, is subtracted from the mea- 
sured value of E^°. The typical size of this correction is 
a few percent of the photon transverse energy. After this 
correction, E^° for truly isolated photons is nominally 
independent of the photon transverse energy. 

In order to make the measurement of E}^° directly com- 
parable to parton-level theoretical predictions, such as 
those described in Section IV, E}^° is further corrected 
by subtracting the estimated contributions from the un- 
derlying event and from pileup. This correction is com- 
puted on an event-by-event basis using a method sug- 
gested in Refs. [43] and [44]. Based on the standard seeds 
for jet reconstruction, which are noise-suppressed three- 
dimensional topological clusters [28], and for two differ- 
ent pseudorapidity regions {\r]\ < 1.5 and 1.5 < < 
3.0), a kr jet-finding algorithm [45, 46], implemented in 
Fast Jet [47], is used to reconstruct all jets without any 
explicit transverse momentum threshold. During recon- 
struction, each jet is assigned an area via a Voronoi tessel- 
lation [48] of the 77 — (/) space. According to the algorithm, 
every point within a jet's assigned area is closer to that 
jet than any other jet. The transverse energy density for 
each jet is then computed from the ratio between the jet 
transverse energy and its area. The ambient transverse 
energy density for the event, from pileup and underlying 
event, is taken to be the median jet transverse energy 
density. Finally, this ambient transverse energy density 
is multiplied by the area of the isolation cone to compute 
the correction to E}^°. 

The estimated ambient transverse energy fluctuates 
significantly event-by-event, reflecting the fluctuations in 
the underlying event and pileup activity in the data. The 
mean correction to the calorimeter transverse energy in 
a cone of radius R ~ 0.4 for an event with one pp interac- 
tion is around 440 MeV in events simulated with PYTHIA 
and 550 MeV in HERWIG. In the data, the mean correction 
is 540 MeV for events containing at least one photon can- 
didate with Et > 15 GeV and exactly one reconstructed 
primary vertex, and increases by an average of 170 MeV 
with each additional reconstructed primary vertex. The 
average number of reconstructed primary vertices for the 
sample under study is 1.56. The distribution of measured 
ambient transverse energy densities for photons passing 
the tight selection criteria is shown in Fig. 1. The im- 
pact of this correction on the measured cross section is 
discussed in Section IX B. For a consistent comparison of 
this measurement to a theoretical prediction which incor- 
porates an underlying event model, the method described 
above should be applied to the generated final state in or- 
der to evaluate and apply the appropriate event-by-event 
corrections. 

After the leakage and ambient-transverse-energy cor- 
rections, the E}^° distribution for direct photons in sim- 
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FIG. 1. The measured ambient transverse energy densities, 
using the jet-area method, for events with at least one tight 
photon. The ambient transverse energy contains contribu- 
tions from both the underlying event and pileup. The broad 
distribution reflects the large event-to-event fluctuations. 



FIG. 2. Transverse energy distribution of photon candidates 
selected in data, after reconstruction and preselection (open 
triangles) and after requiring tight identification criteria and 
transverse isolation energy lower than 3 GeV (full circles). 
The photon candidates have pseudorapidity \'rj''\ < f.37 or 
f.52 < |77^| < 1.81. 



VI. SIGNAL EFFICIENCY 



ulated events is centered at zero, with an RMS width of 
around 1.5 GeV (which is dominated by electronic noise 
in the calorimeter). In the following, all photon candi- 
dates with Ei^° < 3 GeV are considered to be experimen- 
tally isolated. This criterion can be related to a cut on 
the transverse isolation energy calculated at the parton 
level in PYTHIA, in order to mimic the isolation criterion 
implemented in JETPHOX. This parton-level isolation is 
the total transverse energy of all partons that lie inside 
a cone of radius R = 0.4 around the photon direction 
and originated from the same quark emitting the photon 
in either ISR or FSR. The efficiency of the experimental 
isolation cut at 3 GeV for photons radiated off partons 
in PYTHIA is close to the efficiency of a parton-level isola- 
tion cut at 4 GeV. This cut on the parton-level isolation 
is equivalent to the same cut on a particle-level isolation, 
which measures the transverse energy in a cone of radius 
R = 0.4 around the photon after hadronisation (with the 
underlying event removed). The experimental isolation 
criterion is expected to reject roughly 50% of background 
candidates with transverse energy greater than 15 GeV. 



About 110 thousand photon candidates satisfy the 
tight identification criteria and have < 3 GeV: 

around 74 thousand are reconstructed as unconverted 
photons and 36 thousand as converted photons. The 
transverse energy distribution of these candidates is 
shown in Fig. 2. For comparison, the initial distribu- 
tion of all photon candidates after the reconstruction and 
preselection is also shown. 



A. Reconstruction and preselection efficiency 

The reconstruction and preselection efficiency, Eroco, is 
computed from simulated events as a function of the true 
photon transverse energy for each pseudorapidity interval 
under study. It is defined as the ratio between the num- 
ber of true prompt photons that are reconstructed - af- 
ter preselection - in a certain pseudorapidity interval and 
have reconstructed E}^° < 3 GeV, and the number of true 
photons with true pseudorapidity in the same pseudora- 
pidity interval and with particle-level transverse isolation 
energy lower than 4 GeV. The efficiency of the require- 
ment > 15 GeV for prompt photons of true transverse 
energy greater than the same threshold is taken into ac- 
count in Section VIII. 

The reconstruction and preselection efficiencies are cal- 
culated using a cross-section-weighted mixture of direct 
photons produced in simulated 7-jet events and of frag- 
mentation photons produced in simulated dijet events. 
The uncertainty on the reconstruction efficiency due to 
the difference between the efficiency for direct and frag- 
mentation photons, and the unknown ratio of the two 
in the final sample of selected signal photons, are taken 
into account as sources of systematic uncertainty in Sec- 
tion IX A. 

The average reconstruction and preselection efficiency 
for isolated prompt photons with Ijy^ruel < 1-37 or 1.52 < 
|777rucl < 1-81 is around 82%; the 18% inefficiency is due 
to the inefficiency of the reconstruction algorithms at low 
photon transverse energy (a few %), to the inefficiency 
of the isolation requirement (5%) and to the acceptance 
loss from a few inoperative optical links of the calorimeter 
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readout [49]. 



B. Identification efficiency 



The photon identification efficiency, eid, is simiiariy 
computed as a function of transverse energy in each, pseu- 
dorapidity region. It is defined as the efficiency for recon- 
structed (true) prompt photons, with measured < 3 
GeV, to pass the tight photon identification criteria de- 
scribed in Section VB. The identification efficiency is 
determined from simulation after shifting the photon 
shower shapes by "shower-shape correction factors" that 
account for the observed average differences between the 
discriminating variables' distributions in data and MC. 
The simulated sample used contains all the main QCD 
signal and background processes. The average differences 
between data and simulation are computed after apply- 
ing the tight identification criteria. The typical size of 
the correction to the MC efficiency is —3% and is always 
between —5% and zero. The photon identification effi- 
ciency after all selection criteria (including isolation) are 
applied is shown in Fig. 3 and in Table I, including the 
systematic uncertainties that are discussed in more detail 
in Section IX A. The efficiencies for converted photons 
are, on average, 3-4% lower than for unconverted photons 
with the same pseudorapidity and transverse energy. 



TABLE I. Isolated prompt photon identification efficiency in 
the intervals of tlie photon pseudorapidity and transverse en- 
ergy under study. 







Identification Efficiency 
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FIG. 3. Efficiency of the tight identification criteria as a func- 
tion of the reconstructed photon transverse energy for prompt 
isolated photons. Systematic uncertainties are included. 



of hits in the silicon trackers. These selection criteria, 
which do not rely on the shape of the electron shower 
in the calorimeter, are sufficient to select a. W ev 
sample with a purity greater than 95%. The identifi- 
cation efficiency of converted photons is taken from the 
efficiency for selected electrons to pass the tight photon 
selection criteria. This approximation is expected to hold 
to within 3% from studies of simulated samples of con- 
verted isolated prompt photons and of isolated electrons 
from W decays. For unconverted photons, the electrons 
in data are used to infer shower-shape corrections. These 
corrections arc then applied to unconverted photons in 
simulation, in order to calculate the unconverted photon 
efficiency from Monte Carlo. The results from the elec- 
tron extrapolation method are consistent with those from 
the simulation, with worse precision due to the limited 
statistics of the selected electron sample. 



As a cross-check, photon identification efficiencies are 
also inferred from the efficiencies of the same identifica- 
tion criteria applied to electrons selected in data from W 
decays. Events containing W —> ev candidates are se- 
lected by requiring: a missing transverse energy greater 
than 25 GeV (corresponding to the undetected neutrino); 
an opening azimuthal angle larger than 2.5 radians be- 
tween the missing transverse energy vector and any en- 
ergetic jets {Et > 15 GeV) in the event; an electron 
transverse isolation energy in a cone of radius 0.4 in the 
r] — 4> space smaller than 0.3 times the electron transverse 
momentum; and a track, associated to the electron, that 
passes track-quality cuts, such as a large amount of tran- 
sition radiation produced in the TRT and the presence 



C. Trigger efficiency 

The efficiency of the calorimeter trigger, relative to the 
photon reconstruction and identification selection, is de- 
fined as the probability for a true prompt photon, passing 
the tight photon identification criteria and with E"^" < 3 
GeV, to pass the trigger selection. It is estimated in two 
steps. First, using a prescaled sample of minimum bias 
triggers, the efficiency of a lower threshold (« 3.5 GeV) 
level- 1 calorimeter trigger is determined. The measured 
efficiency of this trigger is 100% for all photon candidates 
with reconstructed > 15 GeV passing tight identifi- 
cation criteria. Then, the efficiency of the high-level trig- 
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FIG. 4. Photon trigger efficiency, witli respect to recon- 
structed isolated photon passing the tight identification crite- 
ria, as measured in data (circles) and simulated background 
events (triangles). 



ger is measured using the sample of events that pass the 
levcl-1 calorimeter trigger with the 3.5 GeV threshold. 

The trigger efficiency for reconstructed photon candi- 
dates passing tight selection criteria, isolated and with 
E:^ > 15 GeV is found to be = (99.5 ± 0.5)%, con- 
stant within uncertainties over the full Et and 77 ranges 
under study. The quoted uncertainty is obtained from 
the estimation of the possible bias introduced by using 
photon candidates from data, which are a mixture of sig- 
nal and background photon candidates. Using Monte 
Carlo samples the absolute difference of the trigger effi- 
ciency for a pure signal sample and for a pure background 
sample is found to be smaller than 0.5% for isolated tight 



photon candidates with E^ > 15 GeV. 

A comparison between the high-level trigger efficiency 
in data and in the background predicted by the simula- 
tion is shown in Fig. 4. 



Isolation vs. identification sideband counting 
method 



The first technique for measuring the prompt photon 
yield uses the number of photon candidates observed in 
the sidebands of a two-dimensional distribution to esti- 
mate the amount of background in the signal region. The 
two dimensions are defined by the transverse isolation en- 
ergy E^° on one axis, and the photon identification (710) 
of the photon candidate on the other axis. On the iso- 
lation axis, the signal region contains photon candidates 
with E}^° < 3 GeV, while the sideband contains pho- 
ton candidates with E!^° > 5 GeV. On the other axis, 
photon candidates passing the tight identification crite- 
ria ( "tight" candidates) belong to the 71D signal region, 
while those that fail the tight identification criteria but 
pass a background-enriching selection ("non-tight" can- 
didates) belong to the 71D sideband. The non-tight se- 
lection requires photon candidates to fail at least one of 
a subset of the photon tight identification criteria, but to 
pass all criteria not in that subset. All the shower shape 
variables based on the energy measurement in the first 
layer of the electromagnetic calorimeter are used to define 
the background enriching selection, with the exception 
of Ws.tot, since it is found to be significantly correlated 
with the Ei^° of background photon candidates, while 
the photon yield measurement relies on the assumption 
of negligible (or small) correlations between the trans- 
verse isolation energy and the shower shape quantities 
used to define the background enriching selection. 

The signal region (region "A") is therefore defined by 
photon candidates passing the tight photon identifica- 
tion criteria and having experimental E}^" < 3 GeV. The 
three background control regions consist of photon can- 
didates cither: 

• passing the tight photon identification criteria but 
having experimental E}^" > 5 GeV (region "B" ) 



• having E}^° < 3 GeV and passing the background- 
enriching identification criteria (region "C") 



VII. BACKGROUND SUBTRACTION AND 
SIGNAL YIELD DETERMINATION 



A non-negligible residual contribution of background 
candidates is expected in the selected photon sample, 
even after the application of the tight identification and 
isolation requirements. Two methods arc used to esti- 
mate the background contribution from data and to mea- 
sure the prompt photon signal yield. The first one is used 
for the final cross section measurement, while the second 
one is used as a cross check of the former. All estimates 
are made separately for each region of pscudorapidity and 
transverse energy. 



• having E}^" > 5 GeV and passing the background- 
enriching identification criteria (region "D"). 

A sketch of the two dimensional plane and of the sig- 
nal and background control region definitions is shown 
in Fig. 5. 

The method assumes that the signal contamination in 
the three background control regions is small, and that 
the isolation profile in the non-tight regions is the same 
as that of the background in the tight regions. If these 
assumptions hold, then the number of background can- 
didates in the signal region can be calculated by tak- 
ing the ratio of candidates in the two non-tight regions 
{Nc/N]j), and multiplying it by the number of candi- 
dates in the tight, non- isolated region (Nb)- The number 
of isolated prompt photons passing the tight identifica- 
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FIG. 5. Illustration of the two-dimensional plane, defined by 
means of the transverse isolation energy and a subset of the 
photon identification (ID) variables, used for estimating, from 
the observed yields in the three control regions {B, C, D), the 
background yield in the signal region [A). 
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tion criteria is therefore: 

Ni^ = NA-NB^, (1) 

where Na is the observed number of photon candidates 
in the signal region. 

The assumption that the signal contamination in the 
background control regions is small is checked using 
prompt photon MC samples. As the number of signal 
events in the background control regions is always pos- 
itive and non-zero, corrections are applied to limit the 
effects on the final result. For this purpose, Eq. 1 is 
modified in the following way: 

NT = Na~ {Ns ~ -BNf) fj''''''Jl , (2) 
where ck = —1% (for A" G {B,C, D}) are the signal leak- 

^ A 

age fractions extracted from simulation. Typical values 
for cb arc between 3% and 17%, increasing with the pho- 
ton candidate transverse energy; for cc^ between 2% and 
14%, decreasing with E^. is always less than 2%. 
The total effect of these corrections on the measured sig- 
nal photon purities is typically less than 5%. 

The isolated prompt photon fraction measured with 
this method, as a function of the photon reconstructed 
transverse energy, is shown in Fig. 6. The numbers 
of isolated prompt photon candidates measured in each 
pseudorapidity and transverse energy interval are also re- 
ported in Table II. The systematic uncertainties on the 
measured prompt photon yield and fraction in the se- 
lected sample are described in Section IX B. 



FIG. 6. Fraction of isolated prompt photons as a function 
of the photon transverse energy, as obtained with the two- 
dimensional sideband method. 



B. Isolation template fit method 

The second method relies on a binned maximum likeli- 
hood fit to the E}^° distribution of photon candidates se- 
lected in data which pass the tight identification criteria. 
The distribution is fit to the sum of a signal template and 
a background template, determined from control samples 
extracted from data. This is similar to the technique 
employed in [22], but relies less on simulation for signal 
and background templates. The signal template is de- 
termined from the E}^° distribution of electrons from W 
and Z decays, selected using the criteria described in [50]. 
Electrons from W decays are required to fulfill tight selec- 
tion criteria on the shapes of their showers in the electro- 
magnetic calorimeter and to pass track-quality require- 
ments, including the presence of transition-radiation hits. 
They must also be accompanied by E^^^^ > 25 GeV, 
and the electron- i?™"''' system must have a transverse 
mass larger than 40 GeV. Electrons from Z decays are 
selected with looser criteria, but the pair must have an 
invariant mass close to the Z mass. A single signal tem- 
plate is constructed for each region in |?7|; exploiting the 
independence of E}^° from the transverse energy of the 
object (after applying the corrections described in Sec- 
tion VC) to maximize the available statistics. A small 
bias is expected due to differences between the electron 
and photon E^° distributions, especially in regions where 
there is significant material upstream of the calorimeter. 
A shift of the signal template is applied to the electron 
distributions extracted from data to compensate for the 
differences between electrons and photons seen in simu- 
lation. This shift, computed using simulated photon and 
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TABLE II. Observed number of isolated prompt photons in the photon transverse energy and pseudorapidity intervals under 
study. The first uncertainty is statistical, the second is the systematic uncertainty, evaluated as described in Section IX B. 

Isolated prompt photon yield 
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electron samples, mcreases from 100 MeV to 600 MeV 
with increasing \r]'^\. The background template is ex- 
tracted from data for each {Et, bin, using the same 
reverse-cuts procedure as in the two-dimensional side- 
band technique. A simulation-based correction, typically 
of the order of 3-4%, is applied to the final photon frac- 
tion to account for signal which leaks into the background 
template. The fit is performed in each region of \r]'^\ for 
the individual bins in transverse energy, and the signal 
yield and fraction are extracted. An example of such 
a fit is shown in Fig. 7. The results from this alterna- 
tive technique are in good agreement with those from the 
simpler counting method described in the previous sub- 
section, with differences typically smaller than 2% and 
within the systematic uncertainties that are uncorrclated 
between the two methods. 



C. Electron background subtraction 

The background of prompt electrons misidentified as 
photons needs also to be considered. The dominant 
electron production mechanisms are semileptonic hadron 
decays (mostly from hadrons containing heavy flavor 
quarks) and decays of electroweak bosons (the largest 
contribution being from W decays). Electrons from the 
former are often produced in association with jets, and 
have profiles similar to the dominant backgrounds 
from light mesons. They are therefore taken into account 
and subtracted using the two-dimensional sideband tech- 
nique described in Section VII A. Conversely, electrons 
from W and Z decays have profiles that arc simi- 
lar to those of signal photons. The contribution of this 
background to the signal yield computed in Section VII A 
needs therefore to be removed before the final measure- 
ment of the cross section. 

The fraction of electrons reconstructed as photon can- 
didates is estimated from the data, as a function of the 
electron transverse energy and pseudorapidity, using a 
control sample of Z ^ e+e~ decays. The average elec- 
tron misidentification probability is around 8%. Using 



300 



> 
CD 

i 250h 



200- 
150 
100 
50 



• 2010Data,\/s = 7TeV 

Signal Template 
:::] Background Template 
□ Fit Result 

35 < E| < 40 GeV 
jri^l < 0.6 

ATLAS 

j L dt = 880 nb ' 
5 10 15 20 




25 



E!f°[GeV] 

FIG. 7. Example of a fit to extract the fraction of prompt 
photons using the isolation template technique in the region 
< I77I < 0.6 and 35 < < 40 GeV. The signal template 
is derived from electrons selected from W or Z decays, and 
is shown with a dashed line. The background template is de- 
rived from a background-enriched sample, and is represented 
by a dotted line. The estimated photon fraction is 0.85 and 
its statistical uncertainty is 0.01. 



the W ^ ev and Z — > ee cross section times branch- 
ing ratio measured by ATLAS in pp collisions at ^/s = 7 
TeV [50] , the estimated fraction of photon candidates due 
to isolated electrons is found to be on average ~ 0.5%, 
varying significantly with transverse energy. A maximum 
contamination of (2.5% ± 0.8%) is estimated for trans- 
verse energies between 40 and 50 GeV, due to the kine- 
matic distribution of electrons from W and Z decays. 
The uncertainties on these estimates are less than 1% of 
the photon yield. 
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VIII. CROSS SECTION MEASUREMENT 

The differential cross section is measured by comput- 
ing: 

da -/Vyicid U ^g-j 

(/ ^dt) Ai?^ ^trigger Eroco ^ID 

The observed signal yield (A'yicid) is divided by the 
widths of the £'T-intervals (Ai?^) and by the product of 
the photon identification efficiency (eidj determined in 
Section VI B) and of the trigger efficiency relative to pho- 
ton candidates passing the identification criteria (etriggor, 
determined in Section VI C). The spectrum obtained this 
way, which depends on the reconstructed transverse en- 
ergy of the photon candidates, is then corrected for detec- 
tor energy resolution and energy scale effects using bin- 
by-bin correction factors (the "unfolding coefficients" U) 
evaluated using simulated samples. The corrected spec- 
trum, which is then a function of the true photon energy, 
is divided by the photon reconstruction efficiency Ercco 
(Section VIA) and by the integrated luminosity of the 
data sample, / Cdt. 

The unfolding coefficients are evaluated from the ratio 
of the true to reconstructed distributions of photon 
candidates, using PYTHIA isolated prompt photon simu- 
lated samples. This procedure is justified by the small 
bin-to-bin migrations (typically of the order of a few 
%) that are expected, given the good electromagnetic 
calorimeter energy resolution compared to the width of 
the transverse energy intervals used in this analysis (be- 
tween 5 and 40 GeV) . The values of the unfolding coeffi- 
cients are slightly higher than 1 and decrease as a func- 
tion of Et:, approaching 1. They differ from 1 by less 
then 2% in the \rii'^\ region between 0.0 and 0.6, and by 
less than 5-7% in the other two \rp\ regions, where more 
material upstream of the electromagnetic calorimeter is 
present. 



IX. SYSTEMATIC UNCERTAINTIES 

Several sources of systematic uncertainties on the cross 
section are identified and evaluated as described in the 
following sections. The total systematic uncertainty is 
obtained by combining the various contributions, taking 
into account their correlations: uncorrclated uncertain- 
tics arc summed in quadrature while a linear sum of cor- 
related uncertainties is performed. 



A. Reconstruction, identification, trigger 
efficiencies 

The systematic uncertainty on the reconstruction ef- 
ficiency from the experimental isolation requirement is 
evaluated from the prompt photon simulation varying the 



value of the isolation criterion by the average difference 
(of the order of 500 MeV) observed for electrons between 
simulation and data control samples. It is 2.5% in the 
pseudorapidity regions covered by the barrel calorimeter 
and 4.5% in the end-caps. 

The systematic uncertainty on the identification effi- 
ciency due to the photon shower-shape corrections is di- 
vided into two parts. The first term evaluates the impact 
of treating the differences between the distributions of the 
shower shape variables in data and simulation as an av- 
erage shift. This uncertainty is evaluated in the following 
way: 

• A modified description of the detector material 
is used to produce a second sample of simulated 
photon candidates. These candidates have differ- 
ent shower-shape distributions, due to the differ- 
ent amount of material upstream of and within the 
calorimeter. This alternative model contains an ad- 
ditional 10% of material in the inactive volumes 
of the inner detector and 10% of radiation length 
in front of the electromagnetic calorimeter. This 
model is estimated to represent a conservative up- 
per limit of the additional detector material that is 
not accounted for by the nominal simulation. 

• The correction procedure is applied to the nomi- 
nal simulation to estimate the differences between 
the nominal and the alternative simulation. The 
shifts between the discriminating variable distribu- 
tions in the nominal and the alternative simulation 
are evaluated, and are used to correct the shower 
shape variable distributions of the nominal simula- 
tion. 

• The photon efficiency from the nominal simulation 
is recomputed after applying these corrections, and 
compared with the efficiency obtained from the al- 
ternative simulation. 

The difference between the efficiency estimated from the 
nominal simulation (after applying the corrections) and 
the efficiency measured directly in the alternative sample 
(with no corrections) ranges from 3% at E^ ^ 20 GeV 
to less than 1% at E^ ~ 80 GeV. 

The second part of the systematic uncertainty on the 
identification efficiency accounts for the uncertainty on 
the extracted shower-shape correction factors. The cor- 
rection factors were extracted by comparing tight pho- 
tons in data and simulation; to evaluate the uncertainty 
associated with this choice, the same correction factors 
are extracted using loose photons. The difference in the 
final efficiency when applying the tight corrections and 
the loose corrections is then taken as the uncertainty. 
This uncertainty drops from 4% to 1% with increasing 

Additional systematic uncertainties that may affect 
both the reconstruction and the identification efficien- 
cies are evaluated simultaneously for the product of the 
two, to take into account possible correlations. These 
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sources of uncertainty include the amount of material 
upstream of the calorimeter; the impact of pile-up; the 
relative fraction of direct and fragmentation photons in 
data with respect to simulation; the misidentification of a 
converted photon as unconverted; the difference between 
the PYTHIA and HERWIG simulation models; the impact 
of a sporadic faulty calibration of the cell energies in the 
electromagnetic calorimeter; and the imperfect simula- 
tion of acceptance losses due to inoperative readout links 
in the calorimeter. 

Of all the uncertainties which contribute to this mea- 
surement, the largest ones come from the uncertainty on 
the amount of material upstream of the calorimeter (ab- 
solute uncertainties ranging between 1% and 8% and are 
larger at low E^), and from the uncertainty on the identi- 
fication efficiency due to the photon shower-shape correc- 
tions (the absolute uncertainties are in the range 1-5%, 
and are larger at low E^). 

The uncertainty on the trigger efficiency, evaluated as 
described in Section VI C, is 0.5% and is nearly negligible 
compared to all other sources. 



B. Signal yield estimates 

The following sources of systematic uncertainties af- 
fecting the accuracy of the signal yield measurement us- 
ing the two-dimensional sideband technique are consid- 
ered. 



1. Background isolation control region definition 

The signal yield is evaluated after changing the iso- 
lation control region definition. The minimum value of 
E^° required for candidates in the non-isolated control 
regions, which is set to 5 GeV in the nominal measure- 
ment, is changed to 4 and 6 GeV. This check is sensitive 
to uncertainties in the contribution of prompt photons 
from QED radiation from quarks: these photons are less 
isolated than those originating from the hard process. Al- 
ternative measurements are also performed where a max- 
imum value of E}1° is set to 10 or 15 GeV for candidates 
in the non-isolated control regions, in order to reduce the 
correlation between the isolation variable and the shower 
shape distributions seen in simulated events for candi- 
dates belonging to the upper tail of the isolation distribu- 
tion. The largest positive and negative variations of the 
signal yield with respect to the nominal result are taken 
as systematic uncertainties. The signal photon fraction 
changes by at most ±2% in all the transverse energy and 
pseudorapidity intervals. 



2. Background photon identification control region 
definition 

The measurement is repeated reversing the tight iden- 
tification criteria on a number of strip variables rang- 
ing between two (only Fsidc and Wss) and five (all the 
variables based on the first layer of the electromagnetic 
calorimeter). The largest positive and negative variations 
of the signal yield (with respect to the nominal result) 
from these three alternative measurements are taken as 
systematic uncertainties. The effect on the signal pho- 
ton fraction decreases with increasing photon transverse 
energy, and is around 10% for E^ between 15 and 20 
GeV. 



3. Signal leakage into the photon identification background 
control region 

From the photon identification efficiency studies, an 
upper limit of 5% is set on the uncertainty on the fraction 
cc of signal photons passing all the tight photon identi- 
fication criteria except those used to define the photon 
identification control region. The signal yields in each 
I '7''^ I interval are thus measured again after varying 
the estimated signal contamination in the photon iden- 
tification control regions (cp and cd/cb) by this uncer- 
tainty, and the difference with the nominal result is taken 
as a systematic uncertainty. The signal fraction varia- 
tions are always below 6%. 



4. Signal leakage into the isolation background control 
region 

The fractions cb and cjj of signal photons contaminat- 
ing the isolation control regions depend on the relative 
amount of direct and fragmentation photons in the sig- 
nal selected in a certain E'^, \ri'^\ interval, since the latter 
are characterized by larger nearby activity, and therefore 
usually have slightly larger transverse isolation energies. 
In the nominal measurement, the values of cs and cd are 
computed with the relative fractions of direct and frag- 
mentation photons predicted by PYTHIA. A systematic 
uncertainty is assigned by repeating the measurement af- 
ter varying these fractions between 0% and 100%. The 
measured signal photon fraction varies by less than 5%. 



5. Signal photon simulation 

The signal yield is estimated using samples of prompt 
photons simulated with HERWIG instead of PYTHIA to de- 
termine the fraction of signal leaking into the three back- 
ground control regions. The variations of the signal pho- 
ton fractions in each E^, \r]'^\ interval are below 2%. 
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6. Correlations between the isolation and the photon 
identification variables for background candidates 

Non-negligible correlations between the isolation vari- 
able and the photon identification quantities would affect 
Eq. 2: the true number of isolated tight prompt photon 
candidates would be 

N7-N^^I^HNs~csNf)^^^p^£^ (4) 

where i? s ^ — fej- would then be different from 

unity. The simulation of background events shows a 
small but non-negligible correlation between the isolation 
and the discriminating shower shape variables used to de- 
fine the photon identification signal and background con- 
trol regions. The signal yields arc therefore recomputed 
with the formula in Eq. 4, using for the value pre- 
dicted by the PYTHIA background simulation, and com- 
pared with the nominal results. The effect is smaller 
than 0.6% in the \-q^\ < 1.37 intervals and around 3.6% 
for 1.52 < < 1.81. 



7. Transverse isolation energy corrections 

The effects of the correction for the underlying 
event on the estimated signal yield are also investigated. 
As this is an event- by- event correction, it cannot be un- 
folded from the observed cross section. The impact of 
this correction is evaluated by estimating the signal yield, 
with and without the correction applied, for events with 
only one reconstructed primary vertex (to eliminate any 
effects of pileup). The estimated signal yields using the 
uncorrected values of E^°, normalized to the yields de- 
rived using the corrected values, show no trend in 
or 77. Furthermore, the impact on the cross section of 
the event-by-event corrections is equivalent to that of an 
average correction of 540 MeV applied to the transverse 
isolation energies of all photon candidates. Similar stud- 
ies in PYTHIA and HERWIG MC yield identical results. 



C. Unfolding coefficients 

The unfolding coefficients used to correct the measured 
cross section for E^ bin-by-bin migrations are computed 
using simulated samples. There are three sources of un- 
certainties on these coefficients. 



1. Energy scale uncertainty 

The uncertainty on the energy scale was estimated to 
be ±3% in test beam studies [51], and is confirmed to be 
below this value from the comparison of the Z ^ e~^e~ 



invariant mass peak in data and Monte Carlo. The un- 
folding coefficients are thus recomputed using simulated 
signal events where the true photon energy is shifted by 
±3%. The coefficients change by ±10%. This uncer- 
tainty introduces a relative uncertainty of about 10% on 
the measured cross section which is fully correlated be- 
tween the different E^ intervals within each pseudora- 
pidity range. 

2. Energy resolution uncertainty 

The uncertainty on the energy resolution may affect 
bin-by-bin migrations between adjacent i?T bins. Test 
beam studies indicate agreement between the sampling 
term of the resolution between data simulation within 
20% relative. Furthermore, studies of the Z — > e+e~ in- 
variant mass distribution in data indicate that the con- 
stant term of the calorimeter energy resolution is below 
1.5% in the barrel and 3.0% in the end-cap (it is 0.7% 
in the simulation). The unfolding coefficients are thus 
recomputed after the reconstructed energy of simulated 
photons smearing to take into account a 20% relative in- 
crease of the sampling term and a constant term of 1.5% 
in the barrel and 3.0% in the end-cap. The resulting vari- 
ation of the unfolding coefficients is always less than 1%. 
The uncertainty arising from non-gaussian tails of the 
energy resolution function is estimated by recomputing 
the coefficients using a prompt photon simulation where 
a significant amount of material is added to the detector 
model. The variations of the unfolding coefficients are 
smaller than 1% in all the pseudorapidity and transverse 
energy intervals under study. 

3. Simulated photon transverse energy distribution 

The unfolding coefficients, computed in E^ intervals 
of non-ncgligiblc size, depend on the initial E^ distribu- 
tion predicted by PYTHIA. An alternative unfolding tech- 
nique [52] is therefore used, which relies on the repeated 
application of Bayes' theorem to iteratively obtain an 
improved estimate of the unfolded spectrum. This tech- 
nique relies less on the simulated original E't distribu- 
tion of the prompt photons. The differences between the 
cross-sections estimated using the bin-by-bin unfolding 
and the iterative Bayesian unfolding are within 2%, and 
are taken into account as an additional systematic uncer- 
tainty. 

D. Luminosity 

The integrated luminosity is determined for each run 
by measuring interaction rates using several ATLAS sub- 
detectors at small angles to the beam line, with the abso- 
lute calibration obtained from beam position scans [53]. 
The relative systematic uncertainty on the luminosity 
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measurement is estimated to be 11% and translates di- 
rectly into a 11% relative uncertainty on the cross section. 



X. RESULTS AND DISCUSSION 

The measured inclusive isolated prompt photon pro- 
duction cross sections da/dE^ are shown in Fig. 8, 9, 
and 10. They are presented as a function of the photon 
transverse energy, for each of the three considered pseu- 
dorapidity intervals. They arc also presented in tabular 
form in Appendix B. The measurements extend from 
= 15 GcV to E^ = 100 GeV spanning almost three 
orders of magnitude. The data are compared to NLO 
pQCD calculations, obtained with the JETPHOX program 
as described in Section IV. The error bars on the data 
points represent the combination of the statistical and 
systematic uncertainties (summed in quadrature): sys- 
tematic uncertainties dominate over the whole considered 
kinematic range. The contribution from the luminosity 
uncertainty (11%) is shown separately (dotted bands) as 
it represents a possible global offset of all the measure- 
ments. The total systematic uncertainties on the theoret- 
ical predictions arc represented with a solid band. They 
arc obtained by summing in quadrature the contributions 
from the scale uncertainty, the PDF uncertainty (at 68% 
C.L.) and the uncertainty associated with the choice of 
the parton-level isolation criterion. The same quantities 
are also shown, in the bottom panels of Fig. 8, 9, and 
10, after having been normalized to the expected NLO 
pQCD cross sections. 

In general, the theoretical predictions agree with the 
measured cross sections for E^ > 25 GeV. For lower Et 



and in the two pscudorapidity regions \r]''\ < 0.6 and 
0.6 < I??''! < 1.37, the cross section predicted by JETPHOX 
is larger than that measured in data. Such low transverse 
energies at the LHC correspond to extremely small val- 
ues of Xt = 2E^ I -y/fi, where NLO theoretical predictions 
are less accurate. In such a regime the appropriate values 
of the different scales are not clearly defined, and the un- 
certainties associated with these scales in the theoretical 
predictions may not be well modeled by simple varia- 
tions of any one scale about the default value of E'^ [54] . 
As the low-i?,^ region is where the fragmentation com- 
ponent has the most significant contribution to the total 
cross section, the total uncertainty associated with the 
NLO predictions at low E'^ may be underestimated. 



XI. CONCLUSION 
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FIG. 8. Measured (dots) and expected (full line) inclusive 
prompt photon production cross sections, as a function of the 
photon transverse energies above 15 GeV and in the pseudo- 
rapidity range |?7''| < 0.6. The bottom panel shows the ratio 
between the measurement and the theoretical prediction. 
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The inclusive isolated prompt photon production cross 
section in pp collisions at a ccnter-of-mass energy = 7 
TeV has been measured using 880 nb~^ of pp colli- 
sion data collected by the ATLAS detector at the Large 
Hadron Collider. 

The differential cross section has been measured as a 
function of the prompt photon transverse energy between 



FIG. 9. Measured (dots) and expected (full line) inclusive 
prompt photon production cross sections, as a function of the 
photon transverse energies above 15 GeV and in the pseudo- 
rapidity range 0.6 < |ry^| < 1.37. The bottom panel shows 
the ratio between the measurement and the theoretical pre- 
diction. 
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FIG. 10. Measured (dots) and expected (full line) inclusive 
prompt photon production cross sections, as a function of the 
photon transverse energies above 15 GeV and in the pseudo- 
rapidity range 1.52 < \ri''\ < 1.81. The bottom panel shows 
the ratio between the measurement and the theoretical pre- 
diction. 



15 and 100 GeV, in the three pseudorapidity intervals 
\t]''\ < 0.6, 0.6 < \r]^\ < 1.37 and 1.52 < Itj^] < 1.81, esti- 
mating the background from the selected photon sample 
and using the photon identification efficiency measure- 
ment described in this paper. The photon identification 
using the fine granularity of the calorimeters. A photon 
isolation criterion is used, after an in situ subtraction 
of the effects of the underlying event that may also be 
applied to theoretical predictions. 

The observed cross sections rapidly decrease as a func- 
tion of the increasing photon transverse energy, span- 
ning almost three orders of magnitude. The precision of 
the measurement is limited by its systematic uncertainty, 
which receives important contributions from the energy 
scale uncertainty, the luminosity, the photon identifica- 
tion efficiency, and the uncertainty on the residual back- 
ground contamination in the selected photon sample. 

The NLO pQCD predictions agree with the observed 
cross sections for transverse energies greater than 25 
GcV, while for transverse energies below 25 GeV the cross 
sections predicted by JETPHOX are higher than measured. 
However, the precision of this comparison below 25 GeV 
is limited by large systematic uncertainties on the mea- 



surement and on the theoretical predictions at such low 
values of xt = 2E^/,/s. 

The measured prompt photon production cross section 
is more than a factor of thirty larger than that measured 
at the Tevatron, and a factor of lO'' larger than for photo- 
production at HERA, assuming a similar kinematic range 
in transverse energy and pseudorapidity. This will allow 
the extension of the measurement up to energies in the 
TeV range after only a few years of data taking at the 
LHC. 
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Appendix A: Definition of photon identification 
discriminating variables 

In this Appendix, the quantities used in the selection 
of photon candidates, based on the reconstructed energy 
deposits in the ATLAS calorimeters, are summarized. 

• Leakage in the hadronic calorimeter 

The following discriminating variable is defined, 
energy deposited in the hadronic 



based on the 
calorimeter: 



Normalized hadronic leakage 



R 



had 



TTihad 
E'-r 



(Al) 



is the total transverse energy Ei^^'^ deposited 
in the hadronic calorimeter, normalized to the 
total transverse energy of the photon can- 
didate. 

In the |?7| interval between 0.8 and 1.37 the energy 
deposited in the whole hadronic calorimeter is used, 
while in the other pseudorapidity intervals only the 
leakage in the first layer of the hadronic calorimeter 
is used. 

• Variables using the second ("middle") layer 
of the electromagnetic calorimeter 

The discriminating variables based on the energy 
deposited in the second layer of the electromagnetic 
calorimeter are the following: 
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— Middle rj energy ratio 

^7x7 



(A2) 



is the ratio between the sum E^^j of the ener- 
gies of the second layer cells of the electromag- 
netic calorimeter contained in a 3 x 7 rectangle 
in 77 X (measured in cell units), and the sum 
-^7x7 of the energies in a 7x7 rectangle, both 
centered around the cluster seed. 

Middle 4> energy ratio 

TD _ -^3x3 

' ~ Ell, 



(A3) 



is defined similarly to i?,,. behaves very 
differently for unconverted and converted pho- 
tons, since the electrons and positrons gener- 
ated by the latter bend in different directions 
in 4> because of the solenoid magnetic field, 
producing larger showers in the <j) direction 
than the unconverted photons. 

Middle lateral width 



W2 



(A4) 



measures the shower lateral width in the sec- 
ond layer of the electromagnetic calorimeter, 
using all cells in a window 77 x = 3 x 5 mea- 
sured in cell units. 

• Variables using the first ("front") layer of 
the electromagnetic calorimeter 

The discriminating variables based on the energy 
deposited in the first layer of the electromagnetic 
calorimeter are the following: 



Front side energy ratio 

E{±i) - E{±1) 



idc 



£;(±i) 



(A5) 



measures the lateral containment of the 
shower, along the rj direction. E{zLn) is the 
energy in the ±n strip cells around the one 
with the largest energy. 

Front lateral width (3 strips) 



Ws,3 



Y.E^ 



(A6) 



measures the shower width along rj in the first 
layer of the electromagnetic calorimeter, using 
two strip cells around the maximal energy de- 
posit. The index i is the strip identification 
number, «niax identifies the strip cells with the 
greatest energy, Ei is the energy deposit in 
each strip cell. 



Front lateral width (total) 
Ws,tot measures the shower width along 77 in 
the first layer of the electromagnetic calorime- 
ter using all cells in a window A77 x A(/) = 
0.0625 X 0.2, corresponding approximately to 
20 X 2 strip cells in rj x (p, and is computed as 

Ws,3- 

Front second maximum 



AE = [E. 



•SI 

max 



minj 



(A7) 



is the difference between the energy of the 
strip cell with the second greatest energy 
^^"d„iax' ^^'^ t^*^ energy in the strip cell with 
the least energy found between the greatest 
and the second greatest energy i?,fjjj {AE = 
when there is no second maximum). 



— Front maxima relative ratio 



Er 



atio 



l'^^ max 



(A8) 



measures the relative difference between the 
energy of the strip cell with the greatest en- 
ergy -Efstjnax ^^'^ ^^'^ energy in the strip cell 
with second greatest energy E^,}^ (1 when 
there is no second maximum). 



Appendix B: Cross Section Measurements 

Table III, IV and V list the values of the measured iso- 
lated prompt photon production cross sections, for the 
0.00 < \r]^\ < 0.60, 0.60 < \7f\ < 1.37 and 1.52 < |r;T| < 
1.81 regions, respectively. The various systematic un- 
certainties originating from the purity measurement, the 
photon selection and identification efficiency, the photon 
energy scale and the luminosity are shown. The total un- 
certainty includes both the statistical and all systematic 
uncertainties, except for the uncertainty on the luminos- 
ity. 
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TABLE III. The measured isolated prompt photon production cross section, for 0.00 < I?;''! < 0.60. The systematic uncertainties 
originating from the purity measurement, the photon selection, the photon energy scale, the unfolding procedure and the 
luminosity are shown. The total uncertainty includes both the statistical and all systematic uncertainties, except for the 
uncertainty on the luminosity. 

Measured JETPHOX 

jjo- Stat syst syst syst syst syst total da- total 

dE^ dE^ 

^ (purity) (efficiency) (en. scale) (unfolding) (luminosity) uncertainty ^ uncertainty 

[GcV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] 



[15,20) 


5, 


,24 


±0.11 


+0.52 
-0.88 


±0, 


.81 


+0.51 
-0.46 


±0, 


.11 


±0, 


.58 


+ 1.3 

-1.4 


6, 


,8 


+ 1.4 
-0.9 


[20, 25) 


1, 


,88 


±0.05 


+0.18 
-0.20 


±0, 


.21 


+0.14 
-0.14 


±0, 


.04 


±0, 


.21 


±0.36 


2, 


,38 


+0.45 
-0.30 


[25,30) 


0, 


,88 


±0.03 


±0.07 


±0, 


.08 


+0.09 
-0.08 


±0, 


.02 


±0, 


.10 


+0.16 

-0.15 


1, 


,01 


+0.17 
-0.13 


[30,35) 


0, 


,461 


±0.016 


+0.029 
-0.019 


±0, 


.035 


+0.045 
-0.046 


±0, 


.009 


±0, 


.05 


±0.07 


0, 


,50 


+0.10 
-0.04 


[35,40) 


0, 


,254 


±0.011 


+0.017 
-0.015 


±0, 


.019 


+0.027 
-0.025 


±0, 


.005 


±0, 


.028 


±0.04 


0, 


,28 


+0.04 
-0.03 


[40, 50) 


0, 


115 


±0.005 


+0.008 
-0.006 


±0, 


.007 


+0.009 
-0.009 


±0, 


.0023 


±0, 


.013 


+0.017 
-0.016 


0, 


127 


+0.018 
-0.014 


[50, 60) 


0, 


,050 


±0.003 


+0.003 
-0.002 


±0, 


.003 


+0.006 
-0.005 


±0, 


.001 


±0, 


.005 


+0.008 
-0.007 


0, 


,052 


+0.007 
-0.006 


[60,100) 0, 


,0120 


±0.0007 


+0.0007 
-0.0005 


±0, 


.0006 


+0.0013 
-0.0012 


±0, 


.0002 


±0, 


.0013 


+0.0019 
-0.0018 


0, 


,0121 


+0.0014 
-0.0012 



TABLE IV. The measured isolated prompt photon production cross section, for 0.60 < |?7'''| < 1.37. The systematic uncertainties 
originating from the purity measurement, the photon selection, the photon energy scale, the unfolding procedure and the 
luminosity are shown. The total uncertainty includes both the statistical and all systematic uncertainties, except for the 
uncertainty on the luminosity. 

Measured JETPHOX 
da Stat syst syst syst syst syst total da total 

(purity) (efficiency) (en. scale) (unfolding) (luminosity) uncertainty uncertainty 
[GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] [nb/GeV] 



[15,20) 


5.9 


±0.2 


+ 1.8 
-0.5 


±1, 


.0 


+0.6 
-0.5 


±0.1 


±0, 


.6 


+2.3 

-1.4 


8, 


.5 


+ 1.7 
-1.3 


[20, 25) 


2.23 


±0.07 


+0.49 
-0.18 


±0, 


.28 


+0.16 
-0.16 


±0.04 


±0, 


.24 


+0.6 
-0.4 


3, 


.0 


+0.5 
-0.4 


[25,30) 


1.05 


±0.03 


+0.16 
-0.06 


±0, 


.10 


+0.10 
-0.10 


±0.021 


±0, 


.12 


+0.24 
-0.19 


1, 


.28 


+0.18 
-0.16 


[30,35) 


0.52 


±0.02 


+0.06 
-0.03 


±0, 


.04 


+0.05 
-0.05 


±0.011 


±0, 


.06 


+0.11 
-0.09 


0, 


.64 


+0.11 
-0.09 


[35, 40) 


0.313 


±0.014 


+0.029 
-0.021 


±0, 


.024 


+0.035 
-0.032 


±0.006 


±0, 


.034 


+0.06 
-0.05 


0, 


.344 


+0.052 
-0.039 


[40, 50) 


0.146 


±0.006 


+0.014 
-0.011 


±0, 


.009 


+0.013 
-0.013 


±0.003 


±0, 


.016 


+0.025 
-0.022 


0, 


,161 


+0.022 
-0.019 


[50, 60) 


0.062 


±0.004 


+0.005 
-0.004 


±0, 


.003 


+0.006 
-0.006 


±0.001 


±0, 


.007 


+0.010 
-0.009 


0, 


.065 


+0.009 
-0.007 


[60,100) 0.0138 


±0.0008 


+0.0013 
-0.0009 


±0, 


.0007 


+0.0016 
-0.0014 


±0.0003 


±0, 


.0015 


+0.0025 
-0.0022 


0, 


.0154 


+0.0019 
-0.0015 
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TABLE V. The measured isolated prompt photon production cross section, for 1.52 < \rf'\ < 1.81. The systematic uncertainties 
originating from the purity measurement, the photon selection, the photon energy scale, the unfolding procedure and the 
luminosity are shown. The total uncertainty includes both the statistical and all systematic uncertainties, except for the 
uncertainty on the luminosity. 













Measured 






JETPHOX 




da 


Stat 


syst 


syst 


syst 


syst 


syst 


total 


dcr 

dEj, 


total 








(purity) 


(efficiency) (en. scale) (unfolding) (luminosity) 


uncertainty 




uncertainty 


[GeV] 


[nb/GeV] 


[nb/GcV] 


[nb/GcV] 


[nb/GeV] 


[nb/GeV] 


[nb/GeV] 


[nb/GeV] 


[nb/GeV] 


[nb/GeV] 


[nb/GeV] 


[15,20) 


2.9 


±0.1 


+0.8 
-0.3 


±0.5 


+0.3 
-0.3 


±0.1 


±0.3 


+ 1.1 

-0.7 


3.1 


+0.6 
-0.5 


[20, 25) 


1.12 


±0.04 


+0.15 
-0.08 


±0.16 


+0.08 
-0.08 


±0.02 


±0.12 


+0.27 
-0.24 


1.10 


+0.20 
-0.15 


[25, 30) 


0.47 


±0.02 


+0.06 
-0.04 


±0.05 


+0.05 
-0.04 


±0.01 


±0.05 


+0.11 

-0.09 


0.46 


+0.07 
-0.06 


[30,35) 


0.240 


±0.013 


+0.028 
-0.016 


±0.023 


+0.025 
-0.026 


±0.005 


±0.026 


+0.052 
-0.045 


0.233 


+0.037 
-0.030 


[35,40) 


0.142 


±0.009 


+0.018 
-0.010 


±0.012 


+0.014 
-0.013 


±0.0032 


±0.016 


+0.030 
-0.026 


0.126 


+0.020 
-0.015 


[40, 50) 


0.062 


±0.004 


+0.005 
-0.004 


±0.005 


+0.006 
-0.006 


±0.0013 


±0.007 


+0.011 
-0.010 


0.058 


+0.008 
-0.007 


[50, 60) 


0.0237 


±0.0025 


+0.0026 
-0.0028 


±0.0019 


+0.0024 
-0.0022 


±0.0005 


±0.003 


±0.005 


0.0243 


+0.0033 
-0.0027 


[60,100) 0.0066 


±0.0005 


+0.0005 
-0.0003 


±0.0005 


+0.0008 
-0.0007 


±0.0002 


±0.0007 


+0.0013 
-0.0012 


0.0057 


+0.0007 
-0.0006 
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